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■ Abstract 

We present a detailed numerical study of a chaotic classical system and 
its quantum counterpart. The system is a special case of a kicked rotor 
, and for certain parameter values possesses cantori dividing chaotic re- 

gions of the classical phase space. We investigate the diffusion of particles 
through a cantorus; classical diffusion is observed but quantum diffusion 
C , is only significant when the classical phase space area escaping through 

the cantorus per kicking period greatly exceeds Planck's constant. A 
quantum analysis confirms that the cantori act as barriers. We numeri- 
cally estimate the classical phase space flux through the cantorus per kick 
and relate this quantity to the behaviour of the quantum system. We 
introduce decoherence via environmental interactions with the quantum 
f*^ . system and observe the subsequent increase in the transport of quantum 

ly-^ ' particles through the boundary. 

o ' 

On 

C\ ■ 1 Introduction 

Quantum chaos is a relatively new field, comprising the study of the quantized 
versions of systems which are classically chaotic. The current consensus is that 
quantum systems do not show sensitive dependence on initial conditions in the 
same way as classical systems. This would appear to disqualify them from being 
described as chaotic. However, classical chaos is also apparent in other aspects 
of a system's evolution, the quantum analogies of which are of great interest. 
In particular, these include Kolmogorov-Arnold-Moser (KAM) tori and cantori, 
the study of which is the focus of this paper. Our group has previously pub- 
lished experimental studies jl], |[ [| of a particular, classically chaotic system 
containing KAM tori and cantori, which is realized experimentally in the quan- 
tum regime, and simulated both classically and quantum mechanically. This 
paper concentrates on the detailed results of our computer simulations. 

KAM tori and cantori in a classical phase space are predicted to influence 
the corresponding quantum system. An unbroken KAM boundary will prohibit 
classical diffusion through it, while tunneling across the barrier is possible in a 
quantum system. When interaction terms in the perturbing Hamiltonian are 
sufficiently large as to break up the boundary and create a cantorus or turnstile, 
classical particles will quickly diffuse through that cantorus but the quantum 
wavefunction will be inhibited j|, ||, [| . A heuristic model proposes that with the 
presence of a perturbing Hamiltonian, quantum diffusion is constrained when 
the classical phase space area escaping through the cantorus each period is ~ ft 
^ 0| . Even though the barrier has been broken, the quantum wave function 
still appears to tunnel through the cantorus. 

The link between the quantum domain and the familiar classical world re- 
mains a hotly debated topic. The quantum-classical correspondence (QCC) 
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principle requires that quantum mechanics contains the classical macroscopic 
limit. A promising approach to the question of quantum-classical correspon- 
dence is the study of decoherence — the analysis of the effect of coupling to 
the environment, which inevitably occurs in a real system, in terms of quantum 
coherence . We introduce increased environmental interactions into our quan- 
tum simulations in order to test the hypothesis that the resulting behaviour will 
more closely resemble the classical system. 

The structure of this paper is as follows. In Section ^| we introduce the 
specifically designed double-kicked rotor system, and present the results of a 
classical analysis. In Section || we show that this classical evolution can be well 
described by a simple random model. In Section || we study the corresponding 
quantum behaviour and analyze the properties of the system using the Floquct 
method. In Section || we introduce decoherence into the quantum system by 
two different methods, and also generate Wigncr functions in order to help 
understand the origin of classical behavior from a quantum system. Finally a 
summary is contained in Section ^. All the parameters used in our model were 
chosen so as to correspond with those used in our experiments (3, pi . 



2 Classical Double-Kicked Rotor 

The original and most commonly studied system in quantum chaos is the 5- 
kicked rotor. The observation of dynamical localization in the atomic optics 
realization of the i5-kicked rotor ^ [HJ |ll) provided an important experimental 
link to the most studied system in Hamiltonian chaos. However, a periodic 
pulsed potential of finite time duration (as used in the <5-kicked rotor experi- 
ments) produces a KAM boundary, which becomes more noticeable for longer 
pulse widths. If one wishes to study diffusion through a cantorus, a train of 
single pulses is not the best system. The classical phase space outside the first 
long-lived cantorus is not strongly chaotic and contains many regular regions 
which will inhibit particle diffusion. Hence we have studied the dynamics due 
to a train of double pulses. This system has been the subject of experimental 
investigation by our group in its atomic optics manifestation Jl|, ||, |j| . Figure [l] 
displays our double pulse train. We can write the dimensionless form of the 
Hamiltonian as 

2 00 

H = V --K^<\> f(t~n) (1) 

n— — oo 

where p is the dimensionless momentum conjugate to the coordinate <f> and 
f(t) specifies the temporal shape of the pulses. K is the dimensionless 'kicking 
strength' which is the single parameter varied in our investigation of the classical 
system. The double pulse train Hamiltonian can be written as 

2 00 

H = ^ — K a m cos(4> — 2-kitlt) (2) 

m— — oo 

where a m — ^jsinc(^g L ) cos( 2 yjj : ) (with the sine function defined as sinc(a;) = 
sin(a:)/a;). Each pulse is of width a/2 and the leading edge separation of the two 
pulses is given by A. The KAM boundaries at p = ±107r and ±307r correspond 
to zero values for and 0,15. 
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As discussed in pj, below a critical kicking strength K = K*, the phase 
space of the classical system contains KAM tori given approximately by the 
lines p = ±107r for our choice of a = A = 0.1. For K > K* these break up to 
become partially penetrable cantori. Figure || shows the phase space for K = 70. 
We observe chaotic seas on either side of the cantorus around p = 10n. Figure || 
shows the phase space for K = 280. We have stronger chaos with little island 
structure remaining. The KAM torus around p = 30tt is still unbroken for this 
kick strength, and we find that this is the case for all kicking strengths studied in 
this paper. The strongly chaotic seas surrounding the p = IOtt cantorus provide 
an ideal phase space structure for studying the transport of particles through 
the barrier. 

For a system kicked by finite length pulses, numerical solutions for the motion 
are considerably more difficult to generate than for (5-kicks. If the pulses are 
rectangular, then the system alternates between the motion of a pendulum, and 
free rotation. The solutions for pendulum motion arc in terms of Jacobi elliptic 
functions and elliptic integrals. In our simulations, the elliptic integrals are 
efficiently evaluated with a specialized algorithm coded in C |Q, which is then 
interfaced with MATLAB. 

The starting point for the calculation of the momentum distributions is an 
initial distribution uniform in cf>, and Gaussian in p. For our chosen a and A 
we have a p = 3.6ir which gives 99.5% of initial conditions inside p — \10ir\. We 
find that the results are not strongly dependent on initial distribution, provided 
this proportion is close to 1. To obtain numerical distributions, we choose 10 5 
initial conditions randomly from this distribution, and propagate them through 
70 kick cycles. Figure ^ shows examples of our simulated distributions. Once the 
KAM boundaries at p = ±107r are broken the classical particles will eventually 
distribute themselves uniformly between the p = ±307r tori. 

3 Diffusion Model 

A very simple model, which mimics some aspects of the classical momentum 
distributions, is obtained by treating the system as consisting of three regions 
where strong, homogeneous diffusion occurs, divided by two permeable barriers. 
The system is enclosed by unbroken barriers. Applying this model to our system, 
we then assume that after a short time, the distributions both inside and outside 
the permeable barriers are each essentially uniform. We also assume that over 
time, the populations can leak across the penetrable barriers. The system can 
then be described by three variables giving the population associated with each 
region. If we choose our initial condition so that the population in each of the 
two 'outer' regions is small, we expect the populations to decay exponentially 
towards equilibrium values where the overall distribution is completely uniform. 

We now take the barriers to be located at p = ±107r. The outer boundaries 
are at ±307r, so the system is divided into three equal parts. Each has a phase 
space area (in dimensionless units) of 40ir 2 . With each kick, some area of phase 
space F is mapped from inside the cantori to above p = 10-7T. This area is 
referred to as the 'phase space flux per kick' through the p — \Qir cantorus. 
Identical amounts of phase space are mapped from above p — \Qtt to inside, 
from inside to below p = — 107r, and from below p = — 107r back inside. This 
follows from the incompressibility of phase space flow in a Hamiltonian system 
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and the symmetry of this system in p. The assumption of separate uniform 
distributions inside each of the three intervals means that the probability of a 
trajectory inside the cantori being mapped outside is 2F/40ir 2 . The probability 
for a trajectory outside the cantori to be mapped inside is F/AQn 2 . Now let 
P(\p\ < 10n,t) be the probability for a particular trajectory to be inside the 
cantori. Take the initial condition to be P(\p\ < 10%, 0) = 1, so that both 
outside regions are empty. The evolution will maintain the symmetry between 
them so that the probability that the trajectory is in a particular outside region 
is P(p > 10n,t) = P{p < 10n,t) = (1 
distribution caused by one kick is 

P(\p\ < 107r,t + l)-P(|p| < 107r,t) 



We can write this as 

P(\p\ < 107r,t + l)-i= 
By induction, using our initial condition P(\p\ < 107r,0) = 1, we have 

= i + V 

3 3 

where a = ln(l - 3F/407T 2 ) « -3F/40n 2 (a must be small for the model to 
apply). The probability for a trajectory to be outside the cantori is given by 

P(\p\ >107r,t) = H(l- e <**). (8) 

The parameter a is a function of kick strength which we do not know, except 
that it must increase with K. We can reverse the problem, estimating the 
flux F by fitting an exponential decay to the proportion of trajectories outside 
p = ±107r versus t for a particular kick strength. This estimate is informative 
for comparing classical and quantum results. Figure Jq(a) shows P(\p\ > IOtt, t) 
versus t, for a range of kicking strengths while Figure ||(b) shows ln[2/3— P(\p\ > 
10ir,t)], which should be a straight line according to the simple diffusion model 
(after a few kicks). 

We see that the model appears to be reasonably good for K < 280. For 
higher values of K it may be that the transport across the barrier is too fast 
for equilibrium to approximately hold in each region. Also, in the real system, 
the equilibrium population outside will differ slightly from 2/3, so that the 
straight line form of our graph will break down for P(\p\ < IOtt) very close to 
equilibrium. Inspection suggests that the real equilibrium value differs from 2/3 
by < e~ 3 w 0.05. Figure || shows estimates of the flux per kick F through the 
cantori based on fitting straight lines to the curves in Figure H. 



P(\p\ < 107r,t))/2. The change in the 



2F 

--^P(\P\ < IOtt, t) + ^(1 - P(\p\ < IOtt,* 



3F 



F 

1, 
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(4) 



3F \ 1 
— (P(b|<107r,t)-- 



40tt 



(5) 



(6) 
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4 Quantum Momentum Distributions and Flo- 
quet States 

For comparison with classical simulations, we use an initial density matrix 



1 / n 2 k 2 



(m\ po |n) = — exp y--^T J <W (9) 

where <j p = 3.6n, A is a normalization constant and k is Planck's constant in 
our dimensionless system. Again the results are not strongly dependent on the 
initial spread. We simulate the system dynamics for up to 70 kicks. Figure [7] 
shows examples of these distributions, for the same kicking strengths as were 
presented in the classical case, with k = 2.6. Figure ||(a) shows the quantum 
probability for the atom to be outside p = \10ir\ versus t, for a range of kicking 
strengths. Figure |(b) shows ln[2/3 - P(\p\ > 10n,t)]. 

Comparing these graphs to Figure [| we see that, for all of these kicking 
strengths, the quantum behaviour deviates qualitatively from the classical after 
less than 20 kicks, and the probability outside the boundary for a quantum 
system always levels off well below the classical equilibrium value. This is con- 
sistent with the premise that the diffusion is suppressed when the phase space 
flux across the cantorus is < h. This KAM localization is distinct from the more 
widely studied dynamical localization El. 

The time evolution operator U is unitary, so that its eigenvalues are of 
the form oij — exp(— iEj/%) where Ej is a (dimensionless) quasi-energy. Its 
eigenstates {quasi-energy states or Floquet states) satisfy 

U\a i ) = e* s *t*\a i ). (10) 

In a basis made up of these states, the evolution operator is diagonal. They 
are therefore equivalent to eigenstates of the Hamiltonian in a time-independent 
system, and if we examine the system only once per kicking cycle, they are effec- 
tively stationary. If we are interested in the state of the system after quantum 
saturation has occurred, we can examine the asymptotic (long-time average) 
momentum distribution, which can be written 



P(n\ Po ) = lim A ^ (n\pt\n) (11) 

t=0 

1 

lim - ^ (n| U'poUV |n) . (12) 



t=0 



By inserting the spectral decomposition of U (representation in terms of the 
Floquet states) we obtain ||] 

P(n\po) = 5Z ( a o\Po \ a 'j) ■ \ ( n \ a j)\ 2 ( 13 ) 

3 

and if our initial condition is a pure state p — \n ) (no| we have 

P(n|n )=53|(n |a 3 )| 2 .|(n|a,)| 2 . (14) 
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This function characterizes the asymptotic mixing between the two momentum 
states. We can generate pseudo-colour plots of P(n\no), to obtain visual repre- 
sentations of the momentum confinement in the quantum system for particular 
choices of K and k. 

Figure || shows a series of these pseudo-colour plots. We have used k = 2.6 
and N = 128, so that % \n ma x\ > 50ir. The plots use a logarithmic colour-scale. 
Black and dark grey regions correspond to vanishingly small probability. Mid 
grey regions indicate a low probability, at the level of quantum tunneling. Pale 
grey and white regions have a sizeable probability density. Inspection of the 
plots shows that their shape appears to be chiefly determined by the classical 
barriers to momentum diffusion. The light line along the line p = po arises from 
each state being mapped onto itself, to some degree. Each plot also shows a 
central light square with \p\ and |po| less than 107r, indicating the strong coupling 
of each state in this range to all the others. Sharp borders on this square indicate 
strong confinement by the cantori. These borders blur as the cantori become 
less effective. For increased kicking strength we eventually see a corresponding 
square associated with the KAM tori at p = ±307r. Figure 0(a) has K = 80 
and we see very strong signatures of the classical barriers. Figure ^(b) and (c) 
show K = 180 and K = 280 respectively. We see significant penetration of the 
p = ±107r cantori, but their effect is still obvious. The p = ±307r KAM tori 
still provide a very strong barrier. Figure ||(d) has K = 400. The effects of 
the p = ±107r cantori have almost disappeared, but penetration through the 
p = ±307r KAM tori is still fairly small. 

Again, we see strong localization which begins to break down for K k 300. 
Referring to Figure [| we see that the classical flux through each cantorus be- 
comes comparable with ~k = 2.6 when K 250. It has previously been ob- 
served jjl3| that this criterion appears to give a reasonable estimate of the kicking 
strength for which strong quantum localization by a cantorus will be destroyed. 



5 Decoherence via Spontaneous Emission 

In order to study decoherence we consider the atomic optics manifestation of 
this double kicked system |jl|, |2[ [|. A periodic potential is created by two 
counter-propagating laser beams, and atoms are subjected to a force due to 
the dipole potential. The dipole potential is derived by neglecting the excited 
state amplitude of the atom jl4|, [ll|. We will now consider the first order 
effects arising from a small non-zero amplitude. Rather than making coherent 
stimulated transitions which is the usual interaction between the atoms and 
the pulsed potential, an atom in the excited state may, with finite probability, 
decay to the ground state by spontaneous emission. We can treat this effect 
as a stochastic process, which is different for each realization. The effect of 
a spontaneous emission event arises from the recoil momentum imparted to 
the atom by both the photon exciting the atom and the photon emitted by 
the atom. The atom recoils with a change in p of uk where — 1 < u < 1 and 
u = j(±l + cos (3), where the upper and lower signs occur with equal probability 
and (3 is the angle which the spontaneously emitted photon makes with the x 
axis. (3 is random, with a distribution which is a sum of dipole distributions 
over the set of possible atomic orientations. To a fairly good approximation 
p0[, u can be treated as uniformly distributed between —1 and 1. We define 
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r\ to be the probability per kicking cycle that the atom undergoes spontaneous 
emission. 

5.1 Density Matrix Calculations 

To include spontaneous emission in our density matrix calculations, we use the 
following approximate technique. For a particular run we choose a probability 
per kicking cycle r\ that a particular atom will spontaneously emit. Equivalcntly 
T) is the proportion of the atoms in the ensemble represented by the density 
matrix which will spontaneously emit in each cycle. We then discretise u, so 
that u = ±1 with equal probability. Therefore the recoil of an atom in state 
|n) will leave it in state \n — 1) or \n + 1). These states are representable by 
the density matrix, unlike those arising from continuous u. Once per kick, the 
following replacement is made: 

(m\ p\n) <— -T) ((m + 1| p \n + 1} + (m — 1| p \n - 1)) + (1 - rf) (m\ p\n) (15) 

where we apply periodic boundary conditions. 

We have also performed Monte Carlo wavefunction simulations which are 
more realistic but considerably more time consuming. In particular they take 
into account the continuous distribution of recoil momenta in the x direction, 
and the fact that spontaneous emission can occur at any time during the laser 
kick. The recoil due to spontaneous emission is continuous between — k and 
k, and alters the 'ladder' on which coherent dynamics occur for the particu- 
lar atom. The approximate way in which we account for spontaneous emission 
in our density matrix calculation produces results which are negligibly differ- 
ent from the Monte Carlo wavefunction calculation, but are significantly more 
computationally efficient. 

5.2 Momentum Distributions 

Figures [l^ and [ll] show momentum distributions versus number of kicks for 
the quantum double-kicked rotor, with spontaneous emission effects included. 
Figure |lO| (a) shows behaviour for K — 80 and 77 = 2%. Comparing this to 
Figure |7|(a), any broadening of the distributions due to the spontaneous emission 
is not obvious. However the 'spiky' nature of the pure quantum version has 
been significantly suppressed. Figure |ll](a) has the same kicking strength and 
r\ = 5%. Here wc still sec little broadening but even stronger suppression of the 
quantum peaks. Figure |l0|(b), with K = 180 and rj = 2%, again shows similarity 
to the pure quantum case (in Figure |7|(b)), except for some suppression of peaks 
and slight broadening. In Figure 0(b) has K — 180 and rj — 5%. We now 
see unmistakable movement of probability into the wings of the distribution, 
qualitatively resembling that in the classical system for K = 180, in Figure |i|(b), 
although not as strong. Figure |l(](c) has K = 280 and 77 = 2%. There is 
much more flow of probability through the cantori than in the pure quantum 
case in Figure 0(b) , although the spikiness of the distribution is still strikingly 
different from the classical case in Figure |(b). For T) = 5% in Figure [ll](c), 
the distributions might be considered to look more classical than quantum, 
although the rate of transport through the boundary does not match the classical 
case. Finally for K — 400, we again see (Figures 0(d) and |ll](d)) a significant 
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increase in the probability flow due to spontaneous emission, accompanied by 
distribution shapes which qualitatively resemble classical behaviour except for 
a somewhat smaller rate of transport. 



5.3 Wigner Functions 

A convenient way to visualize the information represented by the density matrix 
is in the form of a Wigner function. For a discrete, truncated basis we use the 
toroidal Wigner function as defined in n& , 



2JV-1 



w(X k ,Pi,t) = ^2 exp (*'jv 

3=0 



njk\ 1 + /l+j 



I- J 



(16) 



where Pi = (S/2)Z and Xk — irk/N. This gives a Wigner function defined on a 
2N x 2N grid. Averaging over cells of four adjacent points we reduce the grid 
to N x N. This was implemented in Matlab using a fast Fourier transform 
algorithm. 

Figures [l^ through [l9] show a series of three-dimensional graphs of Wigner 
functions. Each function represents the state of the quantum double-kicked sys- 
tem after 20 kicks, with the same initial condition as for our other simulations. 
Figure |l2| has K = 80 and no spontaneous emission. We see that the Wigner 
function is strongly contained by the classical barriers. It has a complicated 
folded shape with some sharp spikes. In Figure ^ we introduce spontaneous 
emission with rate rj = 5%. The Wigner function is qualitatively similar, but 
close comparison shows that it has become somewhat smoother, with the con- 
trast between peaks and troughs being reduced. In Figure [ti] we have K = 180 
and rj = 0. The Wigner function 'spills' over the classical barriers, although the 
total probability outside is small. We chiefly notice that the Wigner function is 
now much more complicated in shape, varying rapidly in position and momen- 
tum. There are significant negative peaks present. The addition of spontaneous 



emission in Figure 15 again serves to suppress this variation. It is not obvious 
from inspection of these graphs, but the cantori localization is also destroyed to 
some degree by finite rj. 

In Figure |l^, with K = 280 and rj = 0, we again see an increase in the com- 
plexity of the Wigner function. At this kicking strength, the effect of the KAM 
tori around p = ±307r has become apparent. The introduction of spontaneous 
emission in Figure [l?] again serves to smooth this variation somewhat, while 
pushing more (quasi) probability into the wings of the function. Finally we ex- 
amine K = 400. For the situation without spontaneous emission (Figure |l^) the 
function is again very complicated and noisy-looking in shape. The main visible 
change with increased kicking strength is the increase in the function near to 
the p = ±307r boundaries. Introduction of spontaneous emission (Figure [l9|) 
again suppresses the rapid variation to some extent. 

In the limit of small k we expect that the Wigner pseudo-probability dis- 
tribution will tend to a classical probability distribution. The Wigner function 
itself cannot be interpreted as a probability distribution because it is not al- 
ways positive. We can argue that the Wigner functions of states of particularly 
quantum character will exhibit this non-positivity strongly. The normalization 
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of our discrete Wigner function is 



JV 



W(ct> k ,pi) = l 



(17) 



fej=i 



but due to non-positivity 



N 



\W(4>k,Pi)\ > l. 



(18) 



We would like to quantify the 'non-classicality' of a given state with a single 
positive number. A possible choice is 



We refer to S as 'quantum strangeness'. In this paper, we will not analyze this 
quantity in any detail. However, for example, a mixed state consisting of two 
Gaussian wave packets centered at (4>,p) = (0, ±16S) has S = 0.1765, while for 
a superposition state with equally weighted components of the two wavepackets 
S = 0.7647. For the initial state we use for our simulations S = 0. In general the 
larger the value of S, the more non-classical the character of the state. Figure po| 
shows this parameter calculated for the state of our system after 20 kicks, for 
several kicking strengths K and spontaneous emission rates rj — 0, 2% and 5%. 

5.4 Decoherence Versus Heating 

An atom which undergoes spontaneous emission receives a random momentum 
kick. This will obviously lead to an increase in the width of the momentum 
distribution, the net effect being some incoherent absorption of energy from 
the laser beams. This effect is referred to as 'heating'. This however is not the 
main physical mechanism behind the increased diffusion we observe due to spon- 
taneous emission. In Figure |2l](a) and |l|(b) we show the additional diffusion 
through p — ±107r over the pure quantum case, for several kicking strengths and 
two values of spontaneous emission rate rj. We note that the contribution of 
heating to the crossing of this boundary will be strongest when the momentum 
distribution is steep near p — ±107r. If the broadening effect was due entirely 
to heating it should be most pronounced for the smaller kicking strengths K 
where the pure quantum distribution is strongly localized by the cantori, and 
for fixed rj, should not increase with increasing K. Figure |2l](a) shows that 
the additional diffusion caused by spontaneous emission rate r\ = 2% is much 
larger at K = 180, 280 and 400, than at K = 80. Figure ||(b) shows the same 
trend for rj = 5%. There must therefore be another mechanism which increases 
the diffusion rate, and which is much stronger than the heating mechanism for 
K > 180. This mechanism is the destruction of quantum coherence, or decoher- 
ence, caused by the randomizing effect of spontaneous emission. We have also 
experimentally verified that heating is negligible for K > 150 fTi| . 



N 




(19) 



fc,i=i 
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5.5 Measurement Decoherence, 
or the Anti-Zeno Effect 

Our model for spontaneous emission leads to decoherence which we can consider 
to be 'environmentally induced'. Kaulakys and Contis |L8| discuss the effect of 
projective momentum measurements on the dynamics of the quantum ^-kicked 
rotor. They find that if a momentum measurement is made after every kick, 
then quantum localization is completely destroyed, and unbounded diffusion 
occurs, with the same diffusion constant as in the classical system. They refer 
to this modification of the dynamics as an anti-Zeno effect. Each measurement 
corresponds to the diagonalization of the density matrix in the momentum rep- 
resentation (i.e. off-diagonal elements become zero), or equivalently, the loss of 
all information about position. This effect is formulated without appealing to a 
collapse of the wave-function, and can in principle be produced experimentally. 

We perform a similar simulation for our system to compare this form of de- 
coherence to the spontaneous emission induced decoherence and to the classical 
motion. Computationally this is very similar to the density-matrix spontaneous 
emission simulations. After each double-kick coherent cycle, the density matrix 
is replaced by a matrix with the same diagonal entries, but all zero off-diagonal 
entries. The cycle then repeats. We note that the Wigner functions for states 
generated by this type of simulation are very simple, being the product of a 
uniform position distribution and the momentum distribution. We therefore do 
not show them, but point out that the 'quantum strangeness' parameter S is 
always zero for these functions. 

5.6 Comparison of Quantum and Classical Results 

To compare our quantum and classical results, we calculate the probability for 
a particle to cross the p — ±10tt cantori, P(\p\ > 107r,i), for classical, pure 
quantum, spontaneous emission and anti-Zeno simulations of the system. Fig- 
ure ^(a) has K = 80. Classical and pure quantum simulations give essentially 
horizontal lines with some fluctuation, because the KAM tori present at this 
kicking strength are effective barriers for both systems. The initial conditions 
determine the level of this line. The quantum simulations with spontaneous 
emission show a gradual leakage and the quantum anti-Zeno simulation gives a 
qualitatively similar result. This suggests that, at this kicking strength, deco- 
herence breaks down the quantum cantorus localization, destroying the corre- 
spondence between the quantum and classical systems. 

In Figure |||(b), we present simulations for K = 180. There is now a clear 
difference between the classical and quantum simulations. In fact they appear to 
differ even for very small t, where the quantum curve initially rises more sharply 
than the classical, before saturation sets in. Adding spontaneous emission de- 
stroys the saturation, and for t > 10 causes the quantum system to qualitatively 
more closely resemble the classical. The anti-Zeno simulation however suggests 
that the limit of large decoherence is again faster transport across the cantori 
than in the classical case. We note that this simulation does not involve heating 
effects, so the transport must be due to decoherence. 

Now shifting our attention to K = 280, in Figure |22|(c), we see a large 
difference between classical and pure quantum behaviour which appears to be 
bridged by the introduction of spontaneous emission. The anti-Zeno calculation 
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now closely corresponds to the classical picture, and we conclude that at this 
kicking strength decoherence does make the quantum system 'more classical'. 

Finally we consider K = 400 in Figure ^(d). The pure quantum case is not 
as far from the classical behaviour as for K — 280, but again decoherence is effec- 
tive in increasing the similarity between the systems. The anti-Zeno calculation 
now gives results that are almost indistinguishable from the classical. The close 
correspondence between the anti-Zeno and classical simulations for this and the 
previous figure can be related to the fact that for these kicking strengths the 
size of the classical flux per kick (in Figure |^) has become comparable to our 
dimensionless Planck's constant k = 2.6. 

Our group has previously published experimental results for K — 280 J2j. 
These show good agreement with the spontaneous emission simulations, espe- 
cially the Monte Carlo runs. With the elimination of some systematic experi- 
mental errors we expect that still better agreement would be achieved. 

Figure ^2] demonstrates that, as we have seen, when rj = 0, quantum satura- 
tion sets in by t — 20 over our entire range of kicking strengths. If we examine 
the r\ = results we see that the quantum strangeness S is very small for 
K = 80, but rapidly increases with increasing kicking strength. In general, we 
expect to see S increase with the chaoticity of the corresponding classical sys- 
tem. Classical systems with strong chaos rapidly develop phase space structure 
on all scales, which cannot be reproduced in a quantum system with finite k 7 so 
the quantum system must begin to exhibit non-classical features. It is interest- 
ing here that when K = 80, S is almost negligible, which seems to correspond 
to the fact that the classical and quantum systems are both strongly localized 
by the KAM tori in the classical system. As K increases, the small scale phase 
space structure making up the cantori must differ in the classical and quantum 
systems, and we see that S rises quickly with K in this regime. Examining 
the effects of introducing spontaneous emission, we see that S is significantly 
reduced even by a rate rj = 2%. Along with the 'more classical' diffusion seen in 
Figure |2^, we can see that the Wigner function also exhibits more 'classicality' 
when we introduce environmental decoherence. 



6 Discussion 

We have analyzed and numerically simulated the classical double-kicked rotor 
system and verified that it possesses KAM tori and cantori which present bar- 
riers to diffusion. We have established the success of a simple diffusion model 
for the system and used it to estimate the phase space flux per kick through 
a cantorus as a function of kicking strength. The double-kicked rotor is an 
interesting chaotic system, which would reward further analysis. One aspect 
would be locating periodic trajectories of the system, especially the series of 
these trajectories converging to the noble KAM torus or cantorus. 

The quantum double-kicked system shows strong localization corresponding 
to classical KAM tori. The system is also localized by classical cantori and does 
not show the sharp transition shown classically. Whereas the classical system 
eventually reaches a uniform distribution in phase space, the quantum system 
saturates with probability still significantly confined by the classical barriers. 
This saturation is confirmed by a Floquet analysis of asymptotic momentum 
distributions. The effect is still obvious even when the size of the classical 
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flux per kick becomes of the order of our dimensionless Planck's constant, k. 
In addition to the saturation, we observe fluctuating peaks in the momentum 
distributions, which contrast strongly to the flat-topped distributions seen clas- 
sically. The quantum transport through the classical boundary more closely 
resembles the classical situation as the kicking strength increases, but examina- 
tion of the Wigner functions shows that the system becomes more non-classical. 
This is consistent with the general theory of quantum chaos which suggests that 
stronger chaos in the classical system accelerates the appearance of quantum co- 
herence effects H, H, ||. 

The modern theory of environmental decoherence states that the interaction 
of quantum systems with the environment is a necessary condition for the ap- 
pearance of classical behaviour in real systems. The traditional semi-classical 
limit, which in our formulation is given by k — > 0, is unsatisfactory because the 
quantum break time becomes infinite only logarithmically, and can be surpris- 
ingly short for a real macroscopic system. After this time arbitrary quantum 
superpositions may arise. 

We have observed the effects of decoherence on our quantum system, using 
two models; dissipation by spontaneous emission caused by the kicking laser 
beam, and the anti-Zeno effect. Where the classical flux per kick through the 
cantorus is finite, decoherence disrupts quantum saturation and the system is 
qualitatively more like the classical version. Examination of the Wigner function 
also indicates that quantum interference effects are suppressed by the decoher- 
ence. The anti-Zeno calculations are a kind of extreme decoherence limit, and 
we have seen that they reproduce the classical dynamics accurately, provided 
that £ is not too large compared to the classical flux per kick. It appears not 
unreasonable to suppose that this correspondence continues indefinitely. 

There are numerous aspects of decoherence in this system which could be 
further investigated. The introduction of noise into the kicking of the system 
is expected to have a similar effect to that of dissipation through spontaneous 
emission |2Cj. Particular states of the system will be 'resistant' to decoherence. 
These are expected to be the 'classical' states in the limit of strong decoherence 
and small k. We would like to determine some of these states and quantify the 
'decoherence times' for other states, for comparison with the quantum break 
time of the coherent system. It would be interesting to quantitatively compare 
the anti-Zeno effect with decoherence via noise and dissipation. 

Our results support the idea that the apparently classical nature of the 
universe arises entirely from quantum mechanics. In the real world the un- 
predictable behaviour of chaotic systems must ultimately arise from quantum 
uncertainty. As further progress is made in the investigation of decoherence, 
there is hope that physics will develop a consistent picture of a smooth transi- 
tion between quantum and classical descriptions of reality. 
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Figure 1: Double pulse showing definitions of a and A. 

Figure 2: Poincare section for double-kicked rotor with K = 70. 

Figure 3: Poincare section for double-kicked rotor with K = 280. The cantori 
at p = ±10n are no longer visible in the phase space at this resolution. 

Figure 4: Classical momentum distribution as a function of number of kicks 
t, for the standard double-kicked system with (a) K = 80, (b)180, (c)280 and 
(d)400. 

Figure 5: (a) Proportion of trajectories outside p = ±107r, versus number of 
kicks t. Kicking strength for each curve is indicated by the legend, (b) Loga- 
rithmic plot of 2/3 — P(\p\ > IOtt), versus number of kicks t. 

Figure 6: Estimate of phase space flux per kick (F) through cantorus versus 
kicking strength K. 

Figure 7: Quantum momentum distribution as a function of number of kicks t, 
with (a) K = 80, (b)180, (c)280 and (d)400. 

Figure 8: (a) Probability outside p — ±107r, versus number of kicks t. Kicking 
strength K for each curve is 120 (dotted), 150 (dashed), 180 (dot-dashed), 210 
(solid), 250 (heavy dotted) and 280 (heavy dashed), (b) Logarithmic plot of 
2/3 — P(\p\ > 107r), versus number of kicks t. Classically these curves are 
almost straight lines. 

Figure 9: Pseudo-colour plot of asymptotic momentum distribution P(p\po) for 
(a) K = 80, (b)180, (c)280 and (d) 400, k = 2.6. Computed using 128 states, 
upon which P(p\p ) is normalized to 1. 

Figure 10: Quantum momentum distribution as a function of number of kicks t 
with i] = 2% (a) K = 80, (b)180, (c)280 and (d)400. 

Figure 1 1 : Quantum momentum distribution as a function of number of kicks t 
with r/ = 5% (a) K = 80, (b)180, (c)280 and (d)400. 

Figure 12: Wigner function after 20 kicks with K = 80 and no spontaneous 
emission. 

Figure 13: Wigner function after 20 kicks with K = 80 and rj = 5%. 
Figure 14: Wigner function after 20 kicks with K — 180 and -q = 0. 
Figure 15: Wigner function after 20 kicks with K = 180 and r\ = 5%. 
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Figure 16: Wigner function after 20 kicks with K — 280 and r\ = 0. 
Figure 17: Wigner function after 20 kicks with K = 280 and r\ = 5%. 

Figure 18: Wigner function after 20 kicks with K = 400 and r\ = 0. 
Figure 19: Wigner function after 20 kicks with K = 400 and r\ = 5%. 

Figure 20: 'Quantum strangeness' S — J2ki=i {\W(<j>k,Pl)\ — W(<f>k,Pi)) versus 
kicking strength K for spontaneous emission rates rj = 0, 2%, 5% as indicated 
in legend. 

Figure 21: (a) Relative effect of spontaneous emission rate 77 = 2% for different 
kicking strengths K = 80 (solid), 180 (dotted), 280 (dashed) and 400 (dot- 
dashed). Proportion of atoms outside p = ±10tt for quantum simulations with 
77 = and 77 = 2% spontaneous emission are calculated and the difference is 
plotted, (b) same as (a) but with 77 = 5%. 

Figure 22: Comparison of classical and quantum simulations with (a) K = 
80, (b) 180, (c) 280 and (d) 400. Probability of finding a given atom outside 
p = ±107r is plotted versus number of kicks t. Coherent quantum evolution 
(light solid line); Quantum evolution with spontaneous emission rates 77 = 2% 
(dashed line) and 5% (dash-dotted line) respectively. Quantum simulation with 
anti-Zeno effect (dotted line). Classical simulation (heavy solid line). 
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